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We consider the simplest extension to the Starobinsky model, by allowing an extra scalar field to 
help drive inflation. We perform our analysis in the Einstein frame and calculate the power spectra 
at the end of inflation to second order in the slow-roll parameters. We find that the model gives 
predictions in great agreement with the current Planck data without the need for fine-tuning. Our 
results encourage current efforts to embed the model in a supergravity setting. 
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I. INTRODUCTION 

The theory of inflation, a time of accelerated expansion of the early universe, was first proposed in the 1980’s and 
since then a myriad of models have been studied to explain the mechanism underlying it. The issue with inflation 
is that it is a phenomenological construct, which needs to be embedded into a fundamental theory. In the very first 
model the period of inflation was driven by quantum corrections to the Einstein-Hilbert Lagrangian [1] . In its simplest 
and most studied version, the Einstein-Hilbert action includes an additional term quadratic in the Ricci-scalar R 
and is usually called Starobinsky-inflation or also R? inflation. Later models of inflation are based on the dynamics 
of scalar fields [2-4]. We refer to [5] and [6] for reviews and list of references on inflationary models. Another way 
of generating inflation is by considering further modihcations of General Relativity. Models include Gauss-Bonnet 
gravity [7-11] and higher-order-polynomial corrections [12] (see also [13] for a recent discussion on higher-order 
corrections in Starobinsky inflation). 

There has been revived interest in Starobinsky inflation after Planck 2013 results [14], which place it in a favourable 
light with respect to cosmological observables. This model of inflation is remarkably consistent with current cosmo¬ 
logical data, in particular with measurements of the anisotropies of the GMB. It predicts a spectral index rig « 0.96 
with little spectral running and a small amount of gravitational waves. 

The idea of generating i?^-inflation from a more fundamental theory has come into focus. Similar effort to what 
was done in Ending a realisation of Higgs inflation in supergravity, either to formulate the theory in the Jordan frame 
(see e.g. [15, 16[) or in the Einstein frame [17], has been invested in Starobinsky inflation. Several papers have tried 
to embed the model into a fundamental framework, such as supergravity (see e.g. [18-26]). Motivation of this work 
was also to extend the model in such a way that a large tensor-to-scalar ratio can be obtained. In the context of 
Supergravity, the inflaton fields are components of chiral multiplets (or are combinations thereof) with the vector 
components are not playing an important role and the universe is assumed to behave like a Friedmann-Robertson- 
Walker (FRW) universe on large scales. Some of the cases considered are multi-field models, in the sense that not 
only one field contributes to the dynamics of inflation. 

The question can be asked, whether Starobinsky inflation preserve its attractive properties in the presence of 
interacting matter fields. There is no reason to assume that gravity is the only driving force behind inflation, so in 
our setup we consider an extension of the Starobinsky model by including a scalar field in the matter sector. These 
have been studied in the past [27] to obtain a period of double inflation. In our paper, we study the system in the 
Einstein frame, calculating the amplitude of the curvature perturbation at the end of inflation, the spectral index and 
the tensor-to-scalar ratio. The Einstein frame analysis allows us to use the formalism presented in [28] to calculate 
the power spectra for a large number of model degrees of freedom and initial conditions without using the full held 
equations. In [28], the power spectra were calculated up to second order in the slow-roll parameter. 

Our motivation is two-fold: firstly, we want to explore the robustness of Starobinsky-type inflation in the presence 
of matter fields, since these fields are not necessarily expected to be dynamically insignificant. If these fields contribute 
to the dynamics of the very early universe, then predictions of the Starobinsky model will potentially be altered and 
our aim is to quantify this further. Secondly, embedding the Starbinsky model in a supergravity framework motivates 
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the existence of more than degree of freedom, which are usually ignored. Thus, we take a phenomenological approach, 
allowing two degrees of freedom to evolve. The theory we consider has two free parameters, as we will discuss below. 

The paper is organised as follows: In the next Section we present the theoretical setup and describe our numerical 
method. In Section III we discuss our findings. We conclude in Section IV. 


II. THEORETICAL SETUP 

The theory we consider is the simplest extension of Starobinsky’s original model of inflation we can imagine. It 
includes an ii^-term in the Einstein-Hilbert action as well as a massive scalar field x in the matter sector. The full 
action reads 


S = 






— + ^R^ 
2k 2 










( 1 ) 


In this equation, k = where Mp^ is the reduced Planck mass. The parameter /i has units [mass]“^. We will 

perform the analysis in the Einstein-frame, which can be obtained by a conformal transformation. First, we rewrite 
the gravitational sector of the action above as 


S'c= 


-(1 + 2 « m -^)- — 


( 2 ) 


The equation of motion for (j) gives (j) = R, thus the two action are equivalent. Now considering the conformal 
transformation 


dfiiy — (3) 

with 

= 1 + (4) 

we obtain from eq. (2) 


with 


Se 
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(g2aV, _ 
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( 5 ) 

( 6 ) 


and a = ^ and where we have defined the canonical normalised field = 1 + 2K/i'0. We now define the mass of 
the scalaron ijj as [27]: 


m 


2 _ 1 
6k/j, 


( 7 ) 


In the following, we will work in natural units, i.e. we set k = 1. This type of theory has been considered in the past 
in [27], working in the Jordan frame. The model discussed in [23], with the choice a; = 1, m = and = 1/12, is 
a special case of our model, where the fields have the same mass. Our model allows the two fields to have different 
masses and m^. 

Following the literature, we define b(^j;) = —2a'tp and derive the equations of motions for the fields in an expanding 
Robertson-Walker space-time (the dot denotes the derivative with respect to cosmic time): 


'll) + 3H'ip + V.^ = b^e^^x^, 

X + {3H + 2b^iP)x + e-^^V^ = 0 , 


( 8 ) 

(9) 
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Einstein’s field equation give 


H 


2Ml L 


;.2 


V’ 


2b ■ 2 

e X 


and 


1 

3M|, 




( 10 ) 

( 11 ) 


To study the perturbations produced during inflation, we will not work with the fields ip and x but perform a 
field-rotation using the degree of freedom along the field trajectory (denoted a) and the degree of freedom orthogonal 
to it (denoted s). The fields are defined by 


with 


da = cosOdip + sinOe^dx 
ds = e’’cos 9dx —sin 9dip, 


sin0 = . 

Y '0^ + e^^x^ 


( 12 ) 

(13) 


(14) 


Cosmological perturbations in this system have been systematically studied in the past, see e.g. [28-30]. In [28], 
a formalism was developed, which allows to calculate the power spectra to second order in the slow-roll parameter. 
Instead of integrating the full perturbation equations, we will use this formalism to calculate the power spectrum. 
The method consists of two steps: Firstly, we evaluate the power spectrum at horizon crossing (the * denotes the time 
of horizon crossing and the slow-roll parameter are defined with respect to a and s): 


V-R* = „ n* (I - 2e* - He* + 4e*7?a-a* + 4e*^i*Sg*ce*)(I + hpT‘^)x 

o7r"^6* 


1 + ^(Se* + 20e* - - 8e*^i*Se*ce * -AQ*)/(a;) 


(15) 


^Q* + Hq* 2e*A 


Q* 


9{x) 


where 


Vs* = „ ,* (1 - 2e* - He* + Ae^ij^a* + 4e*^i*Se*ce*)(I + k‘^T'^)x 

oTT 6^ 

2 

I + -(3e* + 20e* - 8e*77,^„* - 8e*Ci*Sg*ce * -DQ^)f{x) 

,( 2 ^Q* + ^Q* ‘^^*Dq* 

r* 9 3 



(16) 


Aq — SrjcycT — 6e + + lOeryo-o- ~ I8e^ + Ile^iSgCg 

Bq = 5r]as - 3CiSe + - 9e^iSe + - r]^s^iCg + 

Cq = 3r]crs - 3CiSe + Sery^rs - Oe^iSg + rjcraiisl - V<y8pi4 + ii4c9, 
Dq = irjss - 30ice(I + Sg) + Qeijss - 7e^ice(I + si) + 77a<TCice(l + sg) 

+ V<Ts^lS9cl + ^ i(Sg — Cg ). 
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/(cc) = 2 — 7 — ln2 — Ina;, (18) 

6g{x) = 16 + - 447 + 127 ^ + 247 In 2 - 44 In 2 + 12 In^ 2 

+ 12 In^ a; — 44 In a; + 247 In a; + 24 In a; In 2. (19) 

The next step is to evaluate the power spectrum at the end of inflation. Because of the presence of isocurvature 
perturbations, the power spectrum of the curvature perturbation evolves on superhorizon scales. Defining 
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the final power spectra are given by 


Vn{N) 


Vs{N) 


'Pn* 





2r]^sfi-kn) f dN" 

J N, 


( 22 ) 


(23) 

(24) 


where the denotes the value of the power spectra at horizon crossing and j = D — A. 

Finally, to calculate the power spectrum of tensor perturbations Vj- we use the slow-roll approximation (see e.g. 
[31]): 

J>r = i([l-2(7 + ln2-lM^. (25) 


III. NUMERICAL RESULTS 

The model we consider in this paper has two degrees of freedom (i/) and x) with two free parameters (m,/, and rriy-). 
The initial conditions are specified by the initial values of the fields and their derivatives. We start our fields at zero 
velocity and make sure that inflation lasts long enough, so that the fields are in the slow-roll regime at the time the 
observable scales leave the horizon. We address the following questions: firstly, do the initial conditions of the fields 
have an significant influence on the observables and secondly, what are the restrictions on the masses of the fields in 
light of the latest Planck results? 

Using the formalism presented in the last section, we evaluate the power spectrum and calculate the amplitude Ag 
of the power spectra at the pivot point fepWot = 0.05/i“^Mpc“^, the spectral index and the tensor-to-scalar ratio 
r. To test the robustness of the predictions of the model, we run a large number of different initial conditions for the 
fields Ip and x for ffr® mass ratio Rm = m^/m^ between 0.1 and 5. We find that the variation in the spectral index 
and the tensor-to scalar ratio is < 1%. The tensor-to-scalar ratio r is small, ranging from r « 0.035 to r « 0.07, with 
the details depending on rn,/, and m^. The variation in Ag can be as large as 10%, for given combinations of and 
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i?m, however that seems to be the exception rather than the rule. For the majority of regions in parameter space, the 
variation in As is < 1%. The running of the spectral index is of order 10“^ or smaller and compatible with current 
observations. 

In Figure 1 we illustrate the three types of trajectories that the fields can follow, as they are approaching the global 
minimum at the end of inflation and the associated power spectra. The figure shows that varying the mass ratio Rm 
changes the behaviour of the fields. The heavier ip becomes, the more the model behaves like double inflation, where 
the heavy held drives inhation hrst and then the second held triggers a second period of inhation. In the opposite 
case, where x is heavier, it quickly approaches zero, leaving inhation to be driven by the ip held. When calculating the 
power spectra, we compared the second-order formalism presented in the last section to the results using a hrst-order 
formalism and found very good agreement. This conhrms that for the purpose of this paper a full integration of the 
equations is not necessary. 

The strongest constraint on the parameters to,/, and comes from the amplitude As- In hg. 2 we show the 
amplitude As at the pivot scale for diherent mass ratios Rm = rn^/rriy^. As one can see, the observed amplitude of 
the primordial perturbations puts strong constraint on the masses for a given Rm, with narrower ranges for higher 
mass ratios. 

We also illustrate the results for the spectral index and the tensor-to-scalar ratio in hg. 2. For comparison we plot 
the Icr limit on Ug coming from the Planck 2013 [14] and Planck 2015 publications [32]. Interestingly, the predictions 
of the model are fully consistent with the 2015 Icr results, but not with the 2013 Icr results. When considering the 2cr 
limits, the discrepancy disappears. The constraints on Ug rule out mass ratios larger than 2.5 and considering ratios 
below 0.2 restrains us to the regime of single held inhation driven by ip (see hg. 1). Under variations of the mass ratio 
Rm G [0.2, 2.5] the values predicted by this extension of the Starobinski inhation are in great agreement with current 
observational data. 

The tensor-to-scalar ratio appears to increase with the mass ratio. Given the Icr constraints on Ug, the tensor-to- 
scalar ratio is limited to be r < 0.07, but higher values of r are possible if the limits on Ug are relaxed. 
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FIG. 1: Field trajectories and corresponding power spectra for runs with the same initial conditions and mass ratios Rm ~ 0.2, 
Rm = 1.0 and Rm =5.0 respectively. The spectral indices for the curvatnre power spectrum is too large in these plots, the 
trajectories shown are for illustrative purposes only. 




















































7 



logio("^x) 




. • • ratio = 0.2 

•- -• ratio = 0.5 

•- -• ratio = 1 

- •- -• ratio = 1.5 

•- -• ratio = 2 

•- -• ratio = 2.5 
— Pianck 1 sigma iimit 
- - Pianck 1 sigma iimit 2013 




-6.0 -5.5 -5.0 

logi(,(m^) 



logi(,(m^) 


FIG. 2: Amplitude, spectral index and tensor-to-scalar ratio r of the curvature perturbation power spectra at the pivot scale 
as a function of (in Planck units). 





























IV. CONCLUSIONS 


We considered the simplest extension to the Starobinsky model, where the quantum corrections to gravity and a 
scalar field drive inflation. Our analysis was done in the Jordan frame, where the theory looks like two scalar fields, 
one with canonical and the other with non-canonical kinetic terms, to second order in slow-roll parameters. The great 
agreement between the first and second order analysis confirms that a full integration of the equations is not necessary 
for our analysis. 

We show that there are regions of parameter space which yield realistic inflation without the need for fine-tunning 
(in the sense that for a given m-^ one can find an (or /r) in the range considered here). We showed that the 
model preserves its appealing features; the spectral index is within the limits of the latest observational results, the 
amplitude of the perturbations at the chosen pivot scale is sufficiently small and the tensor-to-scalar ratio is < 0.1. 

We embed our model in the framework proposed in [23], but we later show that by allowing the fields to have 
different masses, you do not necessarily violate the latest constraints coming from the Planck satellite. If this model 
is to be embedded in Supergravity, the two fields don’t necessarily need to belong to the same multiplet. 

In future work we will be considering particle production at the end of inflation. We will analyse the model in the 
context of pre- and reheating to verify the feasibility of the simplest extension to Starobinsky inflation. 
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